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Abstract 

We examine the spectrum of 2-colour lattice QCD with 4 continuum flavours at a finite chemical 
potential (/i) for quark-number, on a 12 3 x 24 lattice. First we present evidence that the system 
undergoes a transition to a state with a diquark condensate, which spontaneously breaks quark 
number at /i = m n /2, and that this transition is mean field in nature. We then examine the 3 
states that would be Goldstone bosons at [i = for zero Dirac and Majorana quark masses. The 
predictions of chiral effective Lagrangians give a good description of the behaviour of these masses 
for [i < m-x/2. Except for the heaviest of these states, these predictions diverge from our mea- 
surements, once \x is significantly greater than m n /2. However, the qualitative behaviour of these 
masses, indicates that the physics is very similar to that predicted by these effective Lagrangians, 
and there is some indication that at least part of these discrepancies is due to saturation, a lattice 
artifact. 
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I. INTRODUCTION 



Recently there has been a renewed interest in the properties o 



nuclear matter — hadronic 
Much of this interest 



ELQ. 



matter at finite non-zero baryon number (and isospin) density 
comes from a reevaluation of the old idea that quark pairs might condense giving rise to a 
transition to a colour superconducting state at high baryon-number density jjj Q. These 
newer studies indicate that the energies associated with this transition are much larger than 
the original estimates so that they could have a significant effect on the equation-of-state of 
nuclear matter. Unfortunately, adding a finite chemical potential for quark/baryon number 
to the Euclidean QCD action renders the fermion determinant complex which precludes 
the naive application of standard lattice simulation methods. (Some advances have been 
made allowing studies at small chemical potential and high temperatures, but their ranges 
of applicability are limited |5j.) For this reason people have turned to the study of models 
which have some of the properties of QCD at finite baryon number but have real, positive 
fermion determinants, allowing lattice simulations. 

One such model is 2-colour QCD with fundamental quarks. With SU(2) co i our , the quarks 
and anti-quarks are in the same representation, leading to a fermion determinant which 
remains real and positive in the presence of a quark-number chemical potential, //. It is 
expected that, at \i — m n /2, this model exhibits a phase transition to a state with a 
colourless diquark condensate which breaks quark number, and associated Goldstone bosons. 
Thus this condensed phase is a superfluid rather than a superconductor. Chiral effective 
Lagrangians have been used to p redict the phase structure of this theory (1 R Q !(]] as 



have random-matrix models 

This predicted phase structure includin g th e mean-field nature of the transition has been 
observed in lattice simulations Q,^, 14 ,1l5 . Iltl. I17I IrS^ . However, the only simulations in 
which it was possible to study the masses of potential Goldstone bosons were performed at a 
quark mass so large that it was difficult to measure those excitations which were significantly 
more massive than the pion Q]. It was also unclear at such large quark masses if the 
masses of these pseudo- Goldstone modes were well separated from the rest of the hadron 
spectrum. Since these earlier simulations were performed at a quark mass of m = 0.1, we 
have performed new simulations at the same /?(1.5) and quark mass m = 0.025 which should 
halve the pion mass m r . These simulations were performed on a 12 3 x 24 lattice again with 
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the quark-number symmetry breaking parameter A = 0.1m and A = 0.2m to keep the explicit 
symmetry breaking, which depends on A/m, small. We have also performed simulations on 
a smaller (8 4 ) lattice to give us some indication of the magnitude of finite size effects. 

In these simulations which had a single staggered quark field corresponding to 4 con- 
tinuum quark flavours, the lattice flavour symmetry atm = A = /i = 0is U(2). When 
this breaks spontaneously, U(2) — > U(l) giving rise to 3 Goldstone bosons. We study the 
behaviour of these 3 states for m ^ 0, /i ^ and A << m, and compare with the leading- 
order predictions from chiral effective Lagrangians. For /j < m 7r /2, the agreement with these 
predictions is excellent. As /x is increased beyond m w /2, the mass of the excitation which 
would be a Goldstone boson for A = and that of the pion, lie consistently below the pre- 
dictions from tree-level effective Lagrangians. Part of this difference might be due to higher 
order corrections in the effective Lagrangian calculations, which cannot be characterized by 
a single parameter in addition to those of tree-level chiral perturbation theory. However, 
at least some of this 'discrepancy' appears to be due to the fact that, at high densities, 
the fermionic constituents of these excitations are revealed, indicating that we are beyond 
the reach of chiral perturbation theory. However, the high density behaviour of the lattice 
theory is strongly affected by saturation, a lattice artifact, so it is unclear how much of this 
would survive to the continuum. 

In section 2 we review the expected pattern of symmetry breaking for lattice 2-colour 
QCD with 1 fundamental staggered quark field (4 flavours), and the (pseudo)-Goldstone 
spectrum associated with this breaking. Throughout our analysis we compare the lattice 
results to analytical calculations derived from effective models. We use a non-linear-sigma 
model (chiral perturbation theory) at leading order similar to that described in [7||, and a 
linear sigma model that models some of the higher order corrections to chiral perturbation 
theory. We present the linear sigma model effective Lagrangian we use to fit our diquark 
condensates, and its predictions for the pseudo-Goldstone spectrum in section 3, while, for 
completeness, the non-linear sigma model is presented in Appendix A. Section 4 describes 
our simulations, presents our results for the condensates and the evidence for mean-field 
scaling. In section 5 we present our measurements of the spectrum of this theory, and 
its comparison with the predictions of effective Lagrangian analyses. Our conclusions are 
presented in section 6. 
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II. SYMMETRY BREAKING IN 2-COLOUR LATTICE QCD 



The staggered fermion action for 2-colour lattice QCD with 1-staggered fermion, i.e. 4 
continuum flavours in the fundamental representation of the colour group is: 

Sf=Y.hm^)+rn}x+\\[x T r2X + Xr2X T ]\ (1) 

sites ^ 1 ' 

where p{fi) is the normal staggered covariant finite difference operator with links in the +t 
direction multiplied by e M and those in the —t direction by e _M . What follows summarizes 
the analysis of the symmetries of this theory presented in detail in 

At n — m — A = 0, this action has a U (2) flavour symmetry which breaks spontaneously 
to U(l). If it breaks forming a chiral condensate, (xx)) there will be 3 broken £7(2) generators 
and 3 Goldstone bosons, namely; 

1 => X^X 

°\ X T T~2X- XT~2X T (2) 

0"2 X T T2X + XT2X T - 

If it breaks forming a diquark condensate, \{x Tr 2X + XJixJ")-, there will again be 3 broken 
generators and 3 Goldstone bosons, 

i x T ^x + x^x T 

&2 ==> XX (3) 
^3 X 1 ^2X - XT2X 1 ■ 

If fi 7^ and m = 0, the staggered fermion actionals invariant under U(l) x /7(1), which is 
completely spontaneously broken by the diquark condensate. Therefore, only the 2 diquark 
states in equation 01 remain Goldstone bosons. 

When m^O, the chiral condensate forms and the 3 Goldstone bosons of equation El gain 
equal masses given by PCAC. As \x is increased from zero, the mass of the pion created by 
\e\ remains constant, since it has zero quark number and does not feel the effect of the 
chemical potential. The energy and hence the mass of the diquark state created by x r 2X T is 
increased to m n + 2/x for diquarks propagating in the +t direction. The mass of the forward 
propagating anti-diquark is decreased to — 2\x. 
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This latter mass vanishes at \x = m n /2, and it becomes a true Goldstone boson. This 
heralds the phase transition to a state in which quark-number is spontaneously broken by a 
diquark condensate, \{x Tt 2X + X T 2X T )- The Goldstone boson is created by the orthogonal 
linear combination, \(x T ' T 2X ~ X T 2X T )i which we note is just that Goldstone mode which is 
common to equations El and 

For fi > m 7r /2, it is useful to introduce the concept of a total condensate S c such that: 

(xx) = S c cosa, (4) 

and 

« (X T T2X + XT~2X T ) = S c sin a. (5) 

We then see that the heaviest of our 3 would-be Goldstone bosons will be that created by 
the operator 

^{X T T2X + XT2X T ) cosa - xxsina, (6) 

which has zero vacuum expectation value. For \i just above the transition, this is predom- 
inantly a diquark/anti-diquark state, while for fi large, it is predominantly the scalar a/fo 
meson, as expected since m can be neglected. This state will be heavy since the scalar 
meson is not a Goldstone boson when fi > 0, even when m — 0. Finally, the pseudo-scalar 
pseudo-Goldstone boson is created by the operator 

XeX cos a - i (x T r 2 ex + X^X T ) sin a (7) 

where effective Lagrangians suggest a — a, with a introduced in (j3J |SJ). For fi small, this 
state is predominantly a pion, while for large [i it is predominantly a pseudoscalar diquark. 
Since this pseudoscalar diquark would be a Goldstone boson if m = 0, its mass should 
approach zero for large fi when m can be neglected. 

When A 0, the remaining Goldstone boson becomes massive. For the details of this A 
dependence and for the fi dependence of the heavier, would-be Goldstone modes above the 
transition, we must turn to effective Lagrangians and chiral perturbation theory. 

III. LINEAR SIGMA MODEL EFFECTIVE LAGRANGIAN 

For small pion mass m^, fi and A, we should be able to use chiral perturbation theory to 
parametrize the behaviour of the condensates, quark-number density and pseudo-Goldstone 
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spectrum of this theory, as was done in 0,Q]. A reworking of this analysis for the symmetries 
of the staggered quark action is presented in the appendix. However, we find that m n , even 
for the lowest quark mass m we use, is too large for tree-level chiral perturbation theory to 
give a quantitative description of the physics of this system except a relatively low chemical 
potential. While next-to-leading order chiral perturbation theory calculations have been 
performed jl]], these do not yet include spectrum calculations. Even the calculation of the 
order parameters has not been extended beyond the neighbourhood of the critical point at 
next-to-leading order. 

We therefore introduce an alternative effective Lagrangian which incorporates at tree 
level some of the properties expected from an all-orders chiral perturbation theory calcu- 
lation. First, it should have the same phase structure and critical exponents (mean field) 
as tree-level chiral perturbation theory, and a critical point at n — At A = 0, the 

spectrum of pseudo-Goldstone bosons for ji < m 7r /2 should be that predicted in the previ- 
ous section from fairly general arguments, and for /i > it should have one massless 



Goldstone boson. The magnitude of the total condensate y (xx) 2 + {\{x Tt 2X + X T 2X T )) 2 
is independent of /i in tree-level chiral perturbation theory, but depends on /i at next-to- 
leading order. Lattice results indeed indicate that the total condensate increases when /i 
increases. Therefore the total condensate should be allowed to vary. Such variation can 
be allowed if the magnitude of the condensate becomes a dynamical field. With the chiral 
perturbation theory Lagrangian, this would be a non-perturbative effect, since it involves 
producing a bound state. Modifying our Lagrangian to explicitly incorporate such excita- 
tions requires replacing the chiral Lagrangian which is of the non-linear sigma model class 
by the corresponding linear sigma model effective Lagrangian. Since we do not intend to use 
this Lagrangian beyond tree level, we do not have to face the problems of trying to formulate 
a chiral perturbation theory based on this Lagrangian ^| . The simplest Lagrangian of this 
form is 

C cS = ^TrV.E^Ej - ^oM^ReTrM^ - -CTrE^f + \t (TrE^) 2 . (8) 



We have used the same conventions as in 



V >1 
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1 
-1 



(9) 



/ 



Ma, 



i sin cf> cos (f> 
cos 4> i sin 



where tan0 — X/m. 

The field £z = (v + a)H contains v, the minimum of the free energy of the linear sigma 
model, as well as the radial and transverse fluctuations around that minimum. The field 
a describes the radial fluctuations around that minimum. The three pseudo-Goldstone 
modes are the transverse fluctuations around that minimum. They are the same as in chiral 
perturbation theory and are contained in the field £ given by: 



(10) 



where 



U = exp 



XT 



with n 



Ps Qr + iQi 
Qr - iQi Ps 



(11) 



and 



E = 



^ % sin a cos a ^ 



cos a i sin a 



(12) 



corresponds to the minimum of the free energy. M n is the pion mass in the presence of the 
Majorana quark mass A, 

2 Vm 2 + A 2 2 

M n = m n, (13) 



where m n is the pion mass at A = /i — 0. The Lagrangian is written in such a way that the 
masses of the pseudo-Goldstone modes are given by M n at zero chemical potential. 

Under a local flavor transformation V G U(2), the different fields transform in the fol- 
lowing way 



E 
B„ 



VY,V T 

VB V V^ - -Vd v V*. 



(14) 
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At \x = and A = 0, the minimum of the free energy corresponds to 



" = V^^- (15) 

In general, the minimum of the free energy is given by minimizing 

£ = ^ - ( v 2 - 4fi 2 v 2 sin 2 a - 2M 2 v v cos(a - 0). (16) 

Comparing this to the phenomenological effective potential, 

£ = - R 4 a R 2 b a 2 R 2 sin 2 a — cmR cos a — cXR sin a. (17) 

used in our earlier work and in the next section, we see that these are identical under the 
substitutions v = Jb/2R/2, ( = 4a/b, and £ = 16/6 2 provided the critical value of fi, 
fi c = m n /2. 

The minimization conditions for the £ of equation EH are 

v(2v/j, 2 sm2a-v M 2 sm(a-(f))) = (18) 
VQ-m^ cos(a — 0) + v(( — 2v £ + 4/x sin a) = 0. 

The computation of the spectrum of the linear sigma model is similar to that for chiral 
perturbation theory given in the appendix. For this case the a field mixes with both the 
Q-modes, but not with the Ps mode. The secular equation for the a and the two Q-modes 
is given by the setting the determinant of the matrix 

E 2 — p 2 + ( — 6t> 2 £ + 4/i 2 sin 2 a 4/i 2 sin 2a 4/jE sin a 

4/i 2 sin 2a E 2 — p 2 + ( — 2v 2 £ + 2/i 2 (l + cos 2a) 4/x£cosa (J9) 
4fiE sin a 4fiE cos a E 2 - p 2 + C - 2v 2 £ + Afi 2 

to zero. The dispersion relation for the Ps modes is 

E 2 = p 2 + ^M 2 cos(a - 0) (20) 

In order to get the masses of the different modes, the secular equation and the dispersion 

relation must be solved together with the minimization equations. The secular equation can 

be cast in a form more similar to that presented in by making use of the relation 

Q n v o sin (b 

C _ 2 ^ + V = M 2 - — — 21 

vq sm a 

which is a consequence of equations El With this form one can see explicitly that in the 
ordered phase for A = 0, where = but sin a ^ 0, there is a true Goldstone mode. 

8 



IV. LATTICE SIMULATIONS AND SCALING 



We have simulated lattice 2-colour QCD with 1 staggered quark field (4 continuum 
flavours) at finite quark-number chemical potential fi on 8 4 and 12 3 x 24 lattices at 
/3 = 4/g 2 = 1.5 (close to the f3 c for N t = 4) and quark mass m = 0.025 in lattice units. Sim- 
ulations were performed with the explicit symmetry-breaking parameter, A = 0.0025, 0.005 
(and zero, for small fi). We used the hybrid molecular-dynamics algorithm, performing 
simulations of 2000 molecular-dynamics time units at each /i and A with dt as small as 
0.0016. The chiral and diquark condensates, the quark number density, and the spectrum 
of candidate pseudo-Goldstone bosons were measured. 

The diquark condensate, \{x Tr 2X + X T 2X T ), is plotted for the larger lattice in figure H^,. 
This condensate is seen to be small for /i < m n /2 = 0.19264(7), and rapidly increases close to 
m n /2. In addition, the decrease with decreasing A suggests that the condensate would vanish 
as A — ► for fi < 0.15, while for fi > 0.225, the condensate appears destined to remain finite 
in this limit. These observations strongly suggest that there is a phase transition, somewhere 
in the range 0.15 < fi < 0.225. To quantify this observation, we have fitted the behaviour 
of these condensates to scaling forms suggested from effective chiral Lagrangians. The fits 
to the simplest form which comes from the tree level analysis of the chiral Lagrangian of 
the non-linear sigma model variety described in the appendix are poor - \ 2 /d.o.f = 64 - 
so we turn again to a form based on the tree level analysis of a Lagrangian of the linear 
sigma model class, as described in section 3, which allows the magnitude of the condensate 
to vary. We have had good experience with such a form in the past for quenched theories 
and for QCD at finite isospin chemical potential 0, Q The diquark condensate is fitted 
to the form 



which derives from equation El where R and a are given from the minimization conditions of 
the previous section. For later reference we note that the prediction for the chiral condensate 
is 




(22) 



(XX) = cR cos a 



(23) 



while the quark-number density is 



j 'o = b /i R 2 sin 2 a. 



(24) 
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The constant c is denned in terms of the critical /i, fj, c , through 

c = b jA^ a + bfl 2 (25) 
m v 

which is equivalent to equation EE provided \x c = m n /2. 

Our best fit to the form of equation |21 has a = 0.662(14), b = 0.736(13), m = 0.02528(4), 
fj, c = 0.19299(9) and y 2 /d.o.f. = 5.4, over the range < \i < 0.4. \i c is in good agreement 
with the measured value of m n /2 = 0.19264(7), while m is close to the value 0.025 used in 
our simulations, considering the quality of the fit. Since the \ 2 f° r this fit appears poor, 
we compare the fit with the measurements both on the 12 3 x 24 lattice which we used 
for the fit, and for the 8 4 lattice in tables IIIIII While the difference between the 8 4 and 







A = 0.0025 






8 4 


12 3 x 4 


fit 




< X T T2X > 


< X T T2X > 


< X T T2X > 


0.000 


0.0730(2) 


0.07351(5) 


0.07361 


0.100 


0.0979(3) 


0.1001(1) 


0.09991 


0.150 


0.1604(9) 


0.1747(3) 


0.17458 


0.175 


0.2448(17) 


0.2847(6) 


0.28645 


0.200 


0.3856(28) 


0.4475(13) 


0.44775 


0.225 


0.5109(44) 


0.5606(13) 


0.56438 


0.250 


0.5900(51) 


0.6320(14) 


0.63519 


0.300 


0.6615(44) 


0.7075(16) 


0.71084 


0.400 


0.7257(50) 


0.7758(16) 


0.77845 


0.500 


0.7410(46) 


0.8098(17) 


0.82255 


0.600 


0.7586(53) 


0.8312(18) 


0.86512 


0.800 


0.6590(45) 






0.900 


0.4236(39) 






1.000 


0.0293(1) 







TABLE I: The diquark condensate 2(x Tt 2X + X T 2X T ) as a function of fj,, for A = 0.0025 

12 3 x 24 'data' suggest that the finite size effects in the measurements on the larger lattice 
are small, they are almost certainly comparable with and probably larger than the statistical 
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A = 0.005 






8 4 


12 3 x 4 


fit 


I 1 


< X T T~2X > 


< X Tt 2X > 


< X Tr 2X > 


0.000 


0.1437(3) 


0.1451(1) 


0.14518 


0.100 


0.1902(5) 


0.1936(2) 


0.19316 


0.150 


0.2893(12) 


0.3030(4) 


0.30342 


0.175 


0.3826(19) 


0.4077(8) 


0.40704 


0.200 


0.4883(23) 


0.5191(8) 


0.51529 


0.225 


0.5820(29) 


0.5962(9) 


0.59698 


0.250 


0.6275(33) 


0.6529(10) 


0.65228 


0.300 


0.6963(33) 


0.7193(12) 


0.71715 


0.400 


0.7509(37) 


0.7817(11) 


0.78057 


0.500 


0.7795(38) 


0.8162(12) 


0.82395 


0.600 


0.7969(40) 


0.8354(12) 


0.86628 


0.800 


0.7078(33) 






0.900 


0.4945(27) 






1.000 


0.0584(2) 







TABLE II: The diquark condensate \{x T ' T 2X + X T 2X T ) as a function of \i for A = 0.005 

errors and the discrepancies between the 'data' and the fits given in these tables, over the 
range of the fits. For this reason we consider the fits to the linear sigma model form to 
be acceptable, indicating that the system undergoes a second order transition with mean- 
field critical exponents at // = m n /2. We have also shown the fit to the tree-level chiral 
perturbation theory mentioned above in our figure. While this clearly has a more limited 
range of validity, it does not appear that unreasonable over this range. However, in order 
to obtain this quality of fit, we were forced to use m = 0.02632(2). This is far enough from 
the true mass m = 0.025 that, as we shall see, the prediction for the chiral condensate is 
considerably poorer. 

Since this fit gives not only the /i and A dependence but also predicts the m dependence, 
we have plotted the predictions of this fit for our old 'data' at m = 0.1 in figure^. 
Considering the fact that m — 0.1 is rather large to expect fits aimed at the chiral limit to 
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SU(2) N ( =4 /S=1.5 m=0.025 12 3 x24 lattice 



A 
X 

N 
h 

~x 

V 



A 

N 

V 



1.0 



0.0 - 




o - A=0.0000 
x - A=0.0025 
o - A=0.0050 



0.0 



0.2 



0.4 



0.6 



SU(2) N f =4 0=1.5 m=0.1 12 d x24 lattice 



1.0 



0.8 



0.6 



0.4 



0.2 




e 



e 



A=0.00 
A=0.01 
A=0.02 



0.0 - 

0.00 0.25 0.50 0.75 1.00 1.25 

fj, 



FIG. 1: The diquark condensate as a function of /x on a 12 3 x 24 lattice. The solid lines are fits to 
the form of equation 1221 The dashed lines are fits to the chiral perturbation theory forms in the 
appendix, a) m = 0.025, b) m = 0.1. 
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work well, the prediction is remarkably good below fi = 0.6 where the effects of saturation 
start to be seen. 

On the smaller lattice we see that the diquark condensate has a broad peak near /x = 0.6, 
beyond which it falls, remaining very small above \i ~ 1. Since, as we shall see later, the 
quark- number density approaches 2, the maximum value allowed by fermi statistics, at these 
fi values, we interpret this fall as a saturation effect, a finite lattice spacing artifact. Further 
evidence that this is indeed a lattice artifact is found by comparing simulations at different 
attice spacings. Comparing the results presented here at (3 — 1.5 with those at (3 = 1.85 
l^ where the lattice spacing is about 2/3 that at (3 = 1.5, we notice that the value of /i 
in lattice units where saturation is reached is consistent with being the same in both cases. 
In addition, the /i of the peak in the diquark condensate for (3 = 1.85 (in lattice units) 
is at least as large as at /3 — 1.5. If these were continuum effects the relevant lattice /1 
values for (3 = 1.85 would be smaller by roughly a factor of 2/3, since /1 = fi p hysicai a - Hence 
we conclude that the decrease in \{x r r 2 x + X T 2X T ) at large \i is purely a lattice artifact. 
Tables HJn] indicate that finite size effects increase with /j,, and that their effect is to depress 
the values of the condensate as /1 increases. Thus we should expect that the infinite lattice 
peak will be at a higher /i value. However, it is probable that this value will still lie below the 
saturation /i, so that the falloff will occur over a range of /i, even on an infinite lattice, rather 
than as an abrupt discontinuity at the saturation value of fi. We come to this conclusion 
based on our experience with the quen ched theory, where we have been able to examine the 
finite size effects more thoroughly 20j . 

In figureEK we show the chiral condensate, (xx) ? as a function of fi. As expected it remains 
approximately constant for /i < m^, above which it falls towards zero. The predictions of 
equation |22l using the parameters obtained from fits to the diquark condensates, are plotted 
on this graph. The agreement appears excellent over the scaling window. For comparison 
we include the predictions for m = 0.1 on our old 'data' at that mass. Here the agreement 
is considerably poorer. Presumably some of this is due to m = 0.1 being too high for 
these chiral effective Lagrangians. In addition, the chiral condensate is expected to be more 
sensitive to cutoff effects, such as saturation, at higher quark mass. Note that the prediction 
from the chiral perturbation theory fit mentioned above does relatively poorly, because 
this predicts that the magnitude of the total condensate is independent of fi which is not 
true. Thus by forcing this form to fit the diquark condensate, the prediction for the chiral 
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condensate must fail. 

Finally we present the quark number density jo a s a function of lattice \x in figure 
for our 12 3 x 24 simulations. Since this data does not extend to the saturation region, we 
also plot the 8 4 'data' in figure EJ along with the (3 = 1.85 results for comparison. Our 
corresponding 'data' from our old m = 0.1 runs is plotted, along with the predictions from 
equation |2H in figure Eb- 

First we note that the 12 3 x 24 measurements start to depart from the predictions from 
our scaling form (equation^)) somewhat earlier than do either of the condensates, namely 
for fi > 0.3. For m = 0.025, this does admit a small scaling window. For m = 0.01, since 
0.3 < TOj/2, there is no scaling window. The fact that the relevant variable appears to be \i 
rather than 2ii/m n) suggests that at least some of this departure is a lattice artifact related 
to saturation. 

As is observed in figure EK the quark-number density rises even more rapidly at larger 
H, until it saturates at 2 near [l — 1. Comparing these (3 = 1.5 measurements with those at 
(3 = 1.85 (figure Eb) indicates that saturation occurs at approximately the same /i in lattice 
units, confirming that saturation is a lattice artifact. However, jo would reach saturation 
at the same lattice /i independent of lattice spacing, if j oc fi 3 , at large \i. This is precisely 
the behaviour expected at large fi in the continuum. (This was pointed out by Son and 



Stephanov for QCD at finite isospin density 



22j.) The linear rise in the diquark condensate 



at large fi, which is a property of the linear-sigma-model fitting forms does predict a cubic 
rise in j , but in every case we have considered the measured jo far exceeds our predictions. 
The hint that this might be real, i.e. not completely attributable to saturation, comes from 
the observation that the onset of this rapid rise in j with \i appears to occur at a larger 
value of fi in lattice units at (3 = 1.5 (figures 0b EK) than at (3 = 1.85 (figure Eb)- This 
is what would be expected if it is a real effect rather than a lattice artifact. This ear- 
lier departure from the predictions from our tree-level analysis of our effective Lagrangian 
should not come as too much of a surprise. This analysis parametrizes the departure from 
tree-level chiral perturbation theory by a single new parameter. This parameter is set by 
our fit to the diquark condensate. Our new Lagrangian retains the same relationship be- 
tween the diquark condensate and the density as the chiral perturbation theory Lagrangian. 
What we are seeing is evidence that this relationship breaks down at a /i value considerably 
less than the saturation value. To retain agreement beyond this point would require more 
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FIG. 2: The chiral condensate as a function of /i on a 12 3 x 24 lattice. The solid lines are from 
the fits to the diquark condensate. The dashed lines are the predictions of the fit to the tree-level 
chiral perturbation theory form, a) m = 0.025, b) m = 0.1. 
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FIG. 3: The quark-number density jo as a function of pi at (3 = 1.5, m = 0.025 on a 12 3 x 24 lattice. 
The solid curves are the predictions from equation 1241 The dashed curves are from tree-level chiral 
perturbation theory, a) m = 0.0025 b) m = 0.1. 
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terms/parameters in our effective Lagrangian. 

V. THE PSEUDO-GOLDSTONE SPECTRUM 

As discussed in section 2 and made quantitative in section 3, spontaneous breaking of 
the lattice U(2) flavour symmetry atm = /i = A = should give rise to 3 Goldstone bosons. 
When these parameters are small, but non-zero, these excitations become pseudo-Goldstone 
bosons, gaining masses dependent on the magnitude of these symmetry-breaking parameters. 
We have measured the connected and disconnected contributions to the propagators of all 
local scalar and pseudoscalar diquarks in our 12 3 x 24 runs every 1 molecular-dynamics time 
step. The connected propagators are measured using noisy estimators of a point source on 
each odd respectively even site/colour of 1 time-slice of the lattice. The disconnected prop- 
agators are calculated using 5 sets of noisy sources defined over the whole odd respectively 
even sublattice, and the noise- diagonal terms are discarded. 

The first state considered is the scalar diquark created by applying the operator \{x Tt 2X~ 
X T 2X T ) to the vacuum. For A = this will be a true Goldstone boson in the diquark 
condensed phase. With the small As we consider, it should have a small mass in the broken 
phase. Tree-level analysis of the chiral perturbation theory Lagrangian discussed in the 
appendix predicts that, at finite A its mass should be given by equation IA7I for the state 
labeled Q. Note that at A = 0, this reproduces the predictions of simpler arguments m G = 
m-x — 2/i for /i < m w /2 and zero for \x > m n /2. For the linear sigma model approach of 
section 3 this is replaced by the lowest lying solution of the secular equation obtained from 
equation At A = the 2 forms are identical. 

We fit our measured propagators P G to the form 

P G (t) = A{exp[-m G t] + exp[-m G (N t - *)]} (26) 

giving the results shown in figure Et- The solid curves in this figure are the predictions from 
equation The dashed curves are from equation I A 71 What we see is that equation 
describes the decrease in mass of this would-be Goldstone boson well for /i < m n /2, and 
the dip near the transition value m n /2. As fi is increased much beyond m w /2, the 'data' 
falls below these predictions. These linear sigma model predictions are a considerable im- 
provement over those of tree level chiral perturbation theory. This gives us confidence, that 
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FIG. 5: The Goldstone boson of quark-number violation as a function of [i. a) m=0.025, b) m=0.1. 
The solid curves are the predictions from our fits. The dashed curves are from tree-level chiral 
perturbation theory. 
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even though these fits fail above the transition, this excitation will still have the expected 
behaviour in the limit A — > 0, indicating that there is a phase transition to a state with a 
diquark condensate, and that this state is the massless Goldstone boson associated with the 
spontaneous breaking of quark number. In figure we compare our predictions with our 
earlier measurements at m = 0.1, and find the agreement to be somewhat worse than at 
m = 0.025. 

Such departures from the predictions of the linear-sigma-model form could well indicate 
that this model of the higher order terms in chiral perturbation theory is too naive to 
correctly predict more than the qualitative nature of the pseudo- Goldstone spectrum. The 
worse agreement for m = 0.1, where in terms of the scaling variable x = 2/i/m^ saturation 
occurs much sooner, suggests that the deviation could be largely due to saturation. 

Next we turn to the consideration of the pseudoscalar pseudo- Goldstone boson. As de- 
scribed in section 2 and in particular in equation [3 at A = this is expected to be the 
pion for /i < m n /2. Above this value it mixes with the pseudoscalar diquark. For large fi it 
should be predominately a pseudoscalar diquark. Since the mixing angle a in equation is 
not known, beyond the predictions of effective Lagrangians which suggest that a = a, we 
choose to fit the pion and pseudoscalar diquark propagators separately to a form analogous 
to that which we used for the Goldstone propagator. The graphs of the measured masses are 
presented in figures 16171 We note that these masses have the expected behaviour in that 
they remain flat from \i — to the vicinity of m n /2, and then commence to fall, becoming 
small for large fi where they should equal the masses of the would-be Goldstone boson. As 
expected these 2 estimators for the mass of the pseudoscalar pseudo-Goldstone boson, are 
in good agreement. 

The effective (chiral) Lagrangian analysis presented in the appendix predicts that this 
pseudoscalar boson should have a mass given by the P s state of equation IA71 while our 
linear-sigma model form is given in equation 1201 We have plotted these curves on our 
'data' in figures IHtCK- Although these predictions have qualitatively the same form as the 
'data', clearly there is no quantitative agreement. (Note that the linear sigma model form of 
section 3 gives a slightly better fit than the chiral perturbation theory form in the appendix.) 
The form of the mass of equation [201 suggests that we compare our data with that of our 
earlier paper where m = 0.1. This we do in figures EblZb- Again we find the falloff in 
mass above is much more rapid for m — 0.1 than for m = 0.025. This comparison 
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FIG. 6: The pion mass as functions of [a, on a 12 3 x 24 lattice at a) m = 0.025, b) m = 0.1. The 
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suggests that at least some of this is the effect of saturation, but if it is all a saturation 
effect, the range over which saturation has an effect is large. Any discrepancy between the 
measurements and fits which remains after removal of the effects of saturation should again 
be taken as an indication that our linear sigma model effective Lagrangian is inadequate to 
describe all departures from tree-level chiral perturbation theory. 

We note that the term in each of our effective Lagrangians which relates jo to the diquark 
condensate is the term proportional to fi 2 . As we see in equation IA6I of the appendix and its 
equivalent for the linear sigma model, this term also contributes to the masses of the pseudo- 
Goldstone bosons. Thus deviations of jo from our predictions, should imply differences in the 
pseudo-Goldstone boson masses from our predictions. The contribution to the state which 
is a true Goldstone boson at A = 0, is proportional to sin 2 a, while that for the pseudoscalar 
is proportional to cos 2 a. Since a rises from zero above fi = m n /2, we would expect any 
departure from our predictions to occur earlier for the pseudoscalar than for the would-be 
Goldstone boson, which is precisely what we see. 

The third pseudo-Goldstone boson is also a scalar state. It is the linear combination of 
the scalar diquark state and the flavour-singlet scalar meson given in equation El Here a 
is unambiguously defined by the requirement that this state has zero vacuum-expectation 
value, so we have calculated the propagator for this state. For // < m n /2 and A = this 
state is a pure diquark state whose propagator is identical to that of the Goldstone state. 
Thus we will find both the Goldstone excitation and the state we want in this propagator. 
This will also remain true at finite A. Since for finite A, there is no phase transition in 
going from fi < m n /2 to fi > m^/2, it follows that these 2 states continue to mix above the 
transition, although the mixing becomes small for fi >> m 7r /2. For this reason we fit our 
propagator to the form 

P s (t) = A{exp[-m G t} + exp[-m G (N t - t)}} + B{exp[-m s t] + exp[-m s (N t - t)}} (27) 

where mj > m G is the mass of this scalar state, in addition to the form with A = 0, which 
is appropriate to the \x = propagator where mg = 777.5 an d to \i » m^/2 where the 
lower lying state decouples. For \x small but non-zero, we extract this mass by fitting the 
propagator of the diquark state obtained by applying the operator x r 2X T f° the vacuum, to 
the form 

P qq = A{exp[-m s t] + ex P [-m G (N t -t))}. (28) 
23 



This is obtained from the observation that for a pure diquark state (as opposed to a mixture 
of a diquark and an antidiquark) at A = 0, the effect of the chemical potential is to add 
2\i to the effective mass for propagation forward in time and subtract 2/i from the effective 
mass for propagation backward in time, thus separating the 2 scalar excitations. The masses 
obtained from these fits are plotted in figure |S1 
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FIG. 8: The scalar mass as functions of n, on a 12 3 x 24 lattice at m = 0.025. The curves are the 
scaling predictions mentioned in the text. 

The analysis of section 3, predicts this mass to be the middle mass obtained from solving 
the secular equation EH while the corresponding prediction from chiral perturbation theory 
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is the mass of given in equation I A 71 of the appendix. We show these predictions on the 
'data' of figure |H1 The agreement is quite good, in contrast to the other 2 masses. In the low 
\x regime, the mass increases as m n + 2/x, as expected from general arguments. Just above 
uin the mass shows a dip past which it resumes its increase, eventually becoming linear 
again with the same slope but zero intercept. Remembering that for m = 0.1, these masses 
should be roughly twice those presented here for m = 0.025, it is easy to see why they were 
too difficult to measure with any precision at that higher quark mass. 

Finally we should mention the radial excitation. Even at A = 0, its mass is > 2.76 for all 
fi. Since this is close to the momentum cutoff (tt) on the lattice, it is not even clear if this 
should be considered as a real state. In any case, this mass is too high to be of more than 
passing interest. 

VI. CONCLUSIONS 

We have simulated 2-colour lattice QCD with one staggered fermion field corresponding 
to 4 continuum flavours, in the fundamental representation (doublet) of the colour group 
(577(2)), at a finite chemical potential /i, and quark mass m = 0.025. As in previous 
simulations, we have observed the transition from the normal state to one characterized by 
a diquark condensate which spontaneously breaks quark-number. We have found that the 
equation-of-state, which describes the dependence of the diquark and chiral condensates on 
fi and the explicit symmetry breaking parameter A, is well approximated by the tree- level 
approximation to a chiral Lagrangian in the linear sigma model class. The critical scaling 
implied by this analysis indicates that the transition is second order with mean-field critical 
exponents, as expected from chiral perturbation theory analyses through next-to-leading 
order. The measured critical value of fi is consistent with m n /2 as expected. Applying 
the predictions of these fits to our earlier simulations at m = 0.1 we see evidence that 
our equation-of-state also gives a reasonable description of the mass dependence of these 
condensates. However, it is clear that saturation effects (a lattice artifact) limit the range of 
applicability of this equation for the higher mass, and we suspect that m = 0.1 might well 
be close to the limit of applicability of this approximation to chiral perturbation theory, if 
not to the applicability of chiral perturbation theory itself. 

Our equation-of-state also predicts the /i, A and m dependence of the quark-number 
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density. Here, it appears that the scaling window is somewhat narrower than that for the 
condensates. While the effect of saturation might explain why this density grows faster 
than the predictions for larger /i, we suspect that we are seeing the limitations of using 
tree-level results from an effective Lagrangian with only one more parameter than the lead- 
ing order chiral perturbation theory Lagrangian, to model departures from tree-level chiral 
perturbation theory. Use of next-to-leading order chiral perturbation theory to fit our mea- 
surements would require going beyond what has already been done |8:j, and going beyond 
next-to- leading order would be extremely difficult. In addition, going to higher order in chi- 
ral perturbation theory introduces more parameters and thus reduces its predictive power. 
At high enough \x chiral perturbation theory will break down. The scale at which chiral 
perturbation theory breaks down is given by the pion decay constant F n . From the nu- 
merous higher order calculations in chiral perturbation theory, this breakdown scale can be 
estimated to be of the order of a few times F w . In our case at /i = 0, we know the chiral 
condensate, the quark mass and the pion mass. From the Gell-Mann-Oakes-Renner rela- 
tion, we find that F n ~ 0.25. We therefore expect that chiral perturbation theory will break 
down when \i ~ 1. However, at \i ~ 0.5, higher-order corrections should already amount to 
20 — 30%. We therefore conclude that it is very doubtful that chiral perturbation theory 
can describe our 'data' for the largest chemical potentials used in our study. Note that this 
F n is the F n for chiral perturbation theory with the lattice symmetries (see appendix). The 
continuum F„ is half this value. 

Our large /x data suggests that jo is increasing consistent with the expected rise cubic in 
/i. Since the arguments leading to this prediction suggest a counting of degrees of freedom 
identical to that for free quarks, this suggests that we are beyond the range of chiral pertur- 
bation theory, since accounting for all these degrees of freedom in terms of hadrons requires 
including hadrons other than the pseudo-Goldstone bosons. Unfortunately it is difficult to 
disentangle such behaviour from the effects of saturation. A fi 3 increase in j is just what 
is required to keep the saturation threshold at a constant \i in lattice units as the lattice 
spacing is varied, which is what we observe. In addition, the decrease of the condensate 
close to saturation, suggests that the quarks are acting like free quarks, which makes it even 
harder to distinguish real effects from saturation induced artifacts. 

The main goal of this project was to measure the spectrum of pseudo-Goldstone bosons 
for this theory, to enhance our knowledge of the pattern of symmetry breaking. Spontaneous 
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breaking of the £7(2) symmetry of the staggered lattice implementation of 2-colour QCD with 
1 staggered quark field corresponding to 4 continuum flavours at m — A = /i — 0, should 
give rise to 3 Goldstone bosons. Fixing m = 0.025, we have studied the variation of these 
now pseudo- Goldstone bosons as functions of \x for 2 choices of A. Comparison with our 
earlier simulations at m — 0.1 gave some indication of the m dependence of these spectra. 
We have obtained the predictions from our linear sigma model effective Lagrangian for these 
masses, using the parameters obtained from our fits to the diquark condensate. 

For /i — our measurements confirm that all 3 pseudo- Goldstone bosons are degenerate 
with mass consistent with the expected M n . m n is approximately proportional to \fm as 
predicted by PCAC. The observed small deviations from PCAC suggest that m — 0.1 is 
beyond the range of the leading order prediction. At A = 0, the lowest mass state is the 
diquark state orthogonal to the condensate. The mass of this state is expected to fall 
linearly to zero as /i is increased to m n /2. Above this phase transition it should remain 
zero, becoming the Goldstone boson of spontaneously broken quark number. At small non- 
zero A, our effective Lagrangian analysis predicts its behaviour. What we have observed is 
that these predictions are good up to and through the transition. As /i is increased much 
beyond this value, the measured mass lies consistently below these predictions. Examining 
the corresponding predictions for the higher quark mass m = 0.1, where the deviation is 
more severe, suggests that at least some of this deviation is coming from saturation, a lattice 
artifact seen when \x is large enough that the Fermi surface approaches the lattice cutoff. 

When A = 0, we expect that the pseudoscalar pseudo- Goldstone mass will remain con- 
stant as \i is increased up to m n /2. Above this it is expected to decrease rapidly, approaching 
zero for large Our measured values at m — 0.025 and non-zero A show evidence for such 
behaviour. However, the rate of decrease in this mass above the transition is significantly 
faster than that predicted by the linear sigma model effective Lagrangian fits. At m — 0.1, 
this decrease is even more precipitous. Here again, there is evidence to suggest that this 
more rapid decrease is due to the effects of saturation. The improvement in going to the 
smaller mass, where the onset of saturation occurs at larger x = 2/x/m — n, supports this in- 
terpretation, and suggests that if the mass were decreased, eventually the predictions would 
fit the 'data'. 

We have noted that the term in the effective Lagrangian which controls the behaviour of 
the quark-number density, also contributes to these pseudo-Goldstone masses, and does so in 
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a way which would be expected to make the discrepancies worse for the pseudoscalar state. 
Hence it is reasonable to assume that the addition of one extra parameter in going from the 
chiral perturbation theory Lagrangian to the linear sigma model Lagrangian, is inadequate 
to parameterize all departures from tree-level chiral perturbation theory The behaviour of 
jo at high n can be obtained from analyses other than those of effective Lagrangians. This 
suggests that one should abandon the use of effective chiral Lagrangians designed for small 
ji and small and adopt a different approach for large ji. 

At A = 0, the final pseudo-Goldstone mass should increase linearly with /i up to m n /2, 
above which it should briefly decrease before continuing its rise. We see evidence for this 
behaviour in our measurements for m = 0.025. The predictions from effective Lagrangians 
are in good agreement with the 'data' for this state. 

We note that some of the other hadrons which could be expected to contribute at high /x, 
are those that would have been pseudo-Goldstone bosons were it not for the flavour symmetry 
breaking of the staggered lattice. On the lattice, the symmetry breaking is U(2) — > U(l), 
giving 3 Goldstone bosons. In the continuum the symmetry breaking is SU(8) — > Sp(8), 
which gives 27 Goldstone bosons. 

To summarize, we have found m = 0.025 to be small enough to see evidence for mean- 
field scaling and to study the spectrum of the 3 pseudo-Goldstone excitations. At this quark 
mass, the pseudo-Goldstone boson masses lie well below those expected for other 'hadrons'. 
However, saturation, where all available fermion levels are filled, is still close enough to the 
phase transition to make it difficult to disentangle real physics from this lattice artifact, even 
though we do find an adequate scaling region for simple observables. For m — 0.1 the scaling 
window is too small to obtain quantitative information. It would be useful to obtain the 
full next-to-leading order analysis of the pseudo-Goldstone spectrum, and the expressions 
for the order parameters beyond leading order in a. However, this is beyond the scope of 
this paper. 
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APPENDIX A: CHIRAL PERTURBATION THEORY AND THE GOLDSTONE 
SPECTRUM 

In this appendix, we construct chiral perturbation theory for the symmetry breaking 
pattern of the staggered fermion action JI} at m = A = /i = 0: U(2) — > £7(1), and we 
study the spectrum of the three Goldstone excitations. These Goldstone modes become 
massive upon the introduction of a nonzero quark mass or a nonzero diquark source. They 
dominate the physics at low energy. In this appendix we study the spectrum in Chiral 
Perturbation Theory. This problem is similar to what can be found in the literature. In 
chiral perturbation theory for Nj quarks in the adjoint representation has been constructed. 
In this case the symmetry breaking pattern is given by SU(2Nf) — > SO(2Nf). Notice that 
for any number of flavors the symmetry breaking pattern of the staggered fermion action is 
U(2N f ) -> SOJ2NA. 

Following % |23j, we construct the low-energy effective Lagrangian for the Goldstone 
modes induced by the spontaneous symmetry breaking U{2) — > U{\). We find that the 
effective Lagrangian is given by 

F 2 

£ eff = — TrV.SV.St - F 2 M*ReTrM^E, (Al) 



where F is the pion decay constant, M 2 = -y/m 2 + \ 2 (ijjip) /2F 2 , and (^V)o is t ne quark- 
antiquark condensate at m = A = ji — 0. 

We have used the same conventions as in . These notations were already introduced in 
section 3, we just need to replace E/ by E in the expressions that appear in section 3, with 

E = UYAJ J \ (A2) 

where 

( iU \ rr ( P s Qr + iQi) 

u = exp Wlth n = ' A3 



and 



i sin a cos a 
cos a i sin a 



(A4) 
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corresponds to the minimum of the free energy with a given by 



Ajj? cos a sin a = M 2 sin (a — 



(A5) 



We now turn to the study of the spectrum. All the computations made in [7( can be 
easily implemented in our case, since our effective Lagrangian is very similar to that studied 
in [fl. 

First at /i = 0, we find that the mass of the three pseudo-Goldstone modes is given by 
M n . At ji ^ 0, we can expand the effective Lagrangian (jAljl to second order in the Goldstone 
fields. 

i—i 

Following pj, we find that the term quadratic in the Goldstone fields in the effective 
Lagrangian is given by 
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The Q and modes are mixed. If we call Q and the eigenmodes, we find that the 
dispersion laws for the different Goldstone modes are given by 
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